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ABSTRACT 


A method is presented for computing the aerodynamic 
influence coefficients (AlCs) for a surface traveling at subsonic 
speed through an oscillatory gust field. The method is based on 
one of the fundamental solutions of unsteady flow theory--that of 
Sears for a two-dimensional! airfoil traveling through a harmonic 


gust field in an incompresaible fluid. 


The harmonic gust AICs relate the aerodynamic control 
point forces to the spanwise variation of the gust amplitude 


through the following definition 
vee) pVvw bs [c.] fw iw} 
B g Fr 4 4 8 


The Aerospace IBM 7090 Computer Program No. HMOdéd provides 


the gust AICs in printed and optional punched-card output formats. 


The program capacity is 100 surface strips and 100 values of 


reduced velocity. 
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SYMBOLS 


Iocal semichord 

Reference semichord 

The odorsen function 

Filement of gust AIC matrix 

Distance between forward and aft control points 
Control point force 

Bessel functions of first kind of orde. zero and one, 
Local reduced frequency 

Reference reduced frequency 

Lift 

Moment 

Surface span 

Velocity 

Reference gust amplitude 

Local gust amplitude 

Leading edge coordinate 

Sweep angle of surface quarter-chord line 

Gust wave length 

A.rmospheric density 


Sears function 


respectively 


iv 


SECTION I 


FORMULATION OF PROBLEM 


A. Introduction 


The random response of a vehicle to atmospheric turbulence can be 
analyzed if the frequency response to a harmonic gust is known. The harmonic 
gust forces acting ona high aspect ratio surface at subsonic speeds can be 
found from the incompressible two-dimensional theory of Sears” from which 
the gust aerodynamic influence coefficients (AICs) may be derived. The pres- 
ent study is an extension of the computational aspects of the method previously 


reported. 
B. Sign Convention 


A consistent sign convention is chosen between the force and gust direc- 
tions. If the gust is specified as an upwash, a positive force acts upward; if 
the gust is specified as a downwash, a positive force acts downward. The 


leading edge coordinate x is chosen as positive aft. 
ins Derivation of Equations 


We define a set of AlCs for a harmonic gust that relate the complex 
amplitudes of the control point forces to the spanwise variation of the gust 


velocity: 
F VW bsl[c /w 1 
a} pv WO 24", g a 


This definition is completely general, being equally as applicable to 


a lifting surface theory as to a strip theory. However, in the case 


of a strip theory, the AICs take ona simplified partitioned form that appears 


as follows: 


oO; 9 LO ttt yo 
| 
0 | "ey 3 0 way | 0 
s 3. &@ '*%. taxes | @ 
Ca aes <a ee (2) 
: | Pee 46, F | 
a ae oe oe ee 
01 0 | 0 leee Bn 
| | 


The first null partition is reserved for control points whose aerodynamic 
forces can be neglected (e.g., external stores) or can be found from some 
other theory (e.g., slender body theory). The remaining partitions are of 
size 2x1 since two control points are necessary on each strip if it has two 


flexible degrees of freedom (assuming a rigid chord). 


To derive the strip AlCs, we assume that the gust lift and moment 
about the airfoil quarter-chord are known. The equivalence between the 
given loads and the control point forces is shown in Fig. I. From Fig. | we 


note the equivalence 


F + F = L (3) 


dF, =--M ; (4) 


o>: (5) 


relationships between the lift 


we must have the 
lo this end we select the incompressible 


To continue the derivation, 
and moment and the gust velocity 
solution of Sears (cf... Ref 1) who gives the lift and moment (when corrected 
tor sweep) as 


Le - 27 cos A pV w bay o(k) exp(-ikx/b) (6) 


(7) 


where the Sears function o(k) is related to the Theodorsen function C(k) as 


O(k) - IG (k) - id (k) | C(k) + iJ, (k) (8) 


The first sweep correction, the factor cos A, adjusts the two-dimensional 

lift curve slope; the second sweep correction, the factor exp(-ikx/b), accounts 
for the phase difference in the gust maximum amplitude reaching the leading 
edge of each surface strip. Substituting Eqs. (6) and (7) into Eq. (5) yields 


the force-gust velocity relationship 


F bAy (k) exp(-ikx/b) 


2n cos A pV w, (9) 


which, by comparison with Fq. (1), yields the harmonic gust AlCs for the 


j'th strip 


kK. -ik.x./b. 
(b./b )(ay,/s) > ( )) exp( ik, x; /b;) 


2n cos A (10) 
0 


C 
Bj 


where the local reduced frequency m is based upon the local semichord length. 


k. = wbh./V (lla) 
J J 

: (2wV/X)(b,/V) (11b) 

= enb./n (llc) 


where \ is the gust wave length, and the reference reduced frequency is 


based upon the reference semichord 


k = 2Wb /X (12) 
r r 


Once the partitions have been obtained for each strip, the total matrix is 


assembled as indicated in Eq. (2). 


D. References 


1. Y. C. Fung. An Introduction to the Theory of Aeroelasticity. 
New York: John Wiley and Sons, Inc., 1955, p. 409. 

2. W. P. Rodden, E. F. Farkas, and F. C. Slack. ''Harmonic Gust 
Aerodynamic Influence Coefficients by Incompressible Strip Theory: 
Analytical Development and Procedure for the IBM 7090 Computer.” Norair 
Division, Northrop Corporation, Report NOR-61-59, 14 April 1961. 


SECTION II 


GENERAL DESCRIPTION OF INPUT 


A. Units 


Since all input is geometrical and the gust matrix is dimensionless, the 


units of length are only required to be consistent--feet or inches (or 


centimeters). 
B. Classes of Numerical Data and Limitations 


Bis Example Problem 


As an example we consider a two-strip wing with the following 


geometrical properties: 


cos A = 0.7500 


b. = £6.00 

s 5.00 ft 
Strip No. A y(ft) b( ft) x( ft) 
l 3.00 3.00 0.20 
2 2.00 1.60 0.30 


We seek the gust AICs for the two reduced velocities 1/k = 1.00 and 5.00, 
and for the steady case (input as 1/k = 0}. 


ya Program Restrictions and Options 


a. The maximum number of wing strips per data deck is 100. 

b. The maximum number of reduced velocities per data deck 
is 100. 

c. Any practical number of input data decks may be stacked 


successively and run during one machine pass. 


A. 


SECTION III 


DATA DECK SETUP 


Loading Order 


Input decks punched from keypunch forms are loaded behind column 


binary deck HMO02. ‘The data for each input deck must be in the following 


order: 


B. 


(1) 
(2) 
(3) 
(4) 
(5) 
(6) 
(7) 


Heading card 
NSTRIP, NFREQ, NPUNCH 
cos A, b, s 
r 
Ay; series 
b. series 
x. series 


(] Ik); series 


Input Data Description 


(1) 


(2) 


The heading card is for data identification; any characters 
desired may be used in Columns 2 through 72. Column 1 should 


be blank. 

Control card (FORMAT 1814) 

(a) NSTRIP = number of strips; < 100 

(b) NFREQ = number of reduced velocities; < 100 


(c) NPUNCH = 0 or blank if aerodynamic matrices are to be 


punched in cards. 


NPUNCH # 0 if no punched output is desired. 


~) 


(3) Constant parameters (FORMAT 6E 12.8) 


(a) cos A = cosine of sweep angle of quarter-chord 
(b) b. - reference semichord 
(c) 4 = surface span 


(4) Ay, series (FORMAT 6E12.8) 


Ay): Ay>,  OY(NSTRIP) - strip widths 
(5) b. series (FORMAT 6E12.8) 
bi b, , b ONSTRIP) local semichords 


(6) x, series (FORMAT 6612.8) 


X10 Xr °° * NSTRIP) leading edge coordinates 


(7) (1/k)); series (FORMAT 6E12.8) 


(1/k), (Ql /k )> reduced velocities 


OK INEREQ © 
NOTE: Each new series starts on a new line (card). 


ae Example Ke ypunch Forms 


Keypunch forms for the example problem are shown on the following 
page. Columns 73 through 80 are for data deck sequencing and may be any 


choice of letters and numbers. 


pages. 


1. 


io 


SECTION IV 


PROGRAM OUTPUT 


Printed Output 


All input data. 


For each reduced velocity: 


a. Reduced velocity value 
b. Elements of AIC matrix partition representing each strip. 
= Punched card identification and sequencing. 


The example problem printed output is shown on the following 
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B. Punched Output 


l. A deck of punched cards (output) from this program is suitable as 


an input deck to other programs requiring the use of AICs. 


oi All punched output is sequenced in order on Columns 73 through 


80 starting with HM0O2000. The data appear in the following order: 
a. Card 1 contains (1/k), (FORMAT 6E12.8). 


b. Card 2 contains m, the size (number of control points) of 
the AIC matrix, and n, the number of strips (partitions) (FORMAT 1814). 


e. The AIC matrix punched in column binary form and its TRA 


card make up the remainder of the punched output for (1/k ),- 

3. The order of Statement 2 is repeated for all reduced velocities 
per input deck. 

4. Fach matrix is punched in compact form by columns. Column 1 


begins in Origin | and Column 2 in Location (1 + matrix size). 


5. The matrix is punched inthe order: Column 1 (real), Column | 
(imaginary); Column 2 (real), Column 2 (imaginary); - ; Column m (real), 


Column m (imaginary). 
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SECTION V 


PROCESSING INFORMATION 


Ope ration 
STANDARD FORTRAN ‘AOR. oo system 


Featimated Machine Time 
ap tame im? uanutes 


NSTRIFP number of strips 


FREQ number of .. duced velocities 

n number of sets (decks) of input data 
n 

T 02+001 


(NSTRIP). x (FREQ). 
pol a J 


Machine Compone nts (Jsed 


About 2000 core storages 
Standard FORTRAN input tape (NT1) 
Standard FORTRAN outpui print tape (NT 2) 


Standard FORTRAN output punch tape (NT3) 
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SECTION VI 


PROGRAM NOTES 


A. Subroutines Used 
RDLN: reads and prints title cards 
AEROPS5: punch AIC matrix 
BESSEL function routine 
BJYO: computes Bessel functions, order zero 
BJY1: computes Bessel functions, order one 
BINPU: binary punch routine 

B. Generalized Tapes 


Input, print, and punch tapes are defined as Units 2, 3, and 12, 
respectively; however, these may be altered by placing the desired units 


on symbolic cards HM020028, HM020029, and HM020030. 


SECTION VII 


FLOW DIAGRAM 


comeure 


osc’. .aToer 
aice 
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SECTION VIII 


SYMBOLIC LISTING 


A partial list of the principal FORTRAN symbole used in the program 


may be related to the physical notation as follows: 


FORTRAN Notation Physical Notation 

BR b. 

B(I) b. 

CGIR (I) Real AICs for strip i 

CGI (1) Imaginary AICs for strip i 
COSLAM cosA 

DELY (1) dy. 

EJOX Bessel functions of the first kind 
EYOX Bessel functions of the second kind 
FREQ (J) l/k.. j'th value 

NSTRIP Number of strips 

NFREQ Number of reduced velocities 

S ~ 


The complete symbolic listing is given on the following pages. 


HARMONIC GUST AERODYNAMIC INFLUENCE COEFFICIENTS 7/iLis 
C BY INCOMPRESSIBLE STRIP THEORY. HMOZO002 
C HMO20003 
DIMENSION CGIRI(100,4), foe aia DELY(100), Ct1G9)s 4MO20004 
1 x(100), ®yvexi100) 
COMMON NIT1, NI2, NT3 HMO 20007 
( HMO20008 
1 FORMAT { L814 ) HMU20009 
2 FORMAT tt G6E12.8 } iy = 
& FORMAT {10 46X, ILSHHARMONIC GUST AICS / IHO 39K, I3, 9H 
i 13+ 20M REDUCED FREQUENCIES / IMO 24K, siaitetaies tee LAMBDA 
2 Fl15.8, 7H, BR = EL5.8, 6H, § = F1S.8 HMO20013 
5 FORMAT { [HO 230, LHI TR, THDELY(I) 15h, SHB) Lox, aetan / HMOZ0014 
l (1H 21K, 13, ELT.8- 2E20.8)) HMO2Z0015 
6 FORMAT C1NO GK, 9H L/KTRD = SELT.G, £4 CEH LSK, SELT.B7) MMOZOOL6 
7? FORMAT (CLHL 21%, 29HHARMONIC GUST AICS. 1/K({R) = ELS. | 
i 21M, NUMBER OF STRIPS = ae 7 
9 FORMAT (1HO 28x, STHALL STRIP PARTITIONS 2 ROWS BY 1 COLUMN, COMPHEOIO019 
LILLE XK ELEMENTS / {HO 14K, SHSTRIP 4k, LOHROW | REAL 9X, LOHROW 1 1MHMO200270 
ZAG 9x, LOWROW 2 REAL 9x, LOHROW 2 IMAG ) 4HM0200?1 
10 FORMAT (1H 5K, LIS, *S619.8 7 MMO20022 
Cc 
C NTL = INPUTIAPE HMO20 
C NY¥2 = OUTPUTTAPE HBO20075 
C NTS = PUNCHTAPE HMO 20026 
U HMOZ0027 
NTi=?2 Je 
NT2 = 3 eS 
NT3 = 1? = 
c 
TwOP!l 6.2831853 HMO20032 
IjTwO = 2 HMO2Z20036 
NC AROD=0 
———. = =e 
14 CALL ROUN CONTE, NT2, DD E: 
READINPUTTAPE NT1, 1, NSTRIP, NFREC, NPUNTH 
READINPUTTAPE NT1, 2, COSLAM, BR, § 


if} 
HA t 
i 
iy 


m™ 
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17 


19 


READINPUTTAPE NTL, 25 (DELVE DT)» T=leNSTRIP) 
READINPUTTAPE NI ly 2, (BCL), Tel »NSTRIP) 
READINPUTTAPE NT1+ 2+ (ACL) s Em LeNSTRIP? 
READINPUTTAPE NIL, 2, CPREQLID, Lai, NFREO) 


CONL = TwOPLeCOSL AM 


WRITEQUTPUTTAPE NI2, 4, NSTRIP, NEREQ, COSL AM, 


WRITEOUTPUTTAPE NI2, S, (CE, Detvtld,. BEI). 
WRITEQUTPUTTAPE NIZ, Gy (FREQULDs T= eNPREQ) 


DO 21 JFREQ = 1, NFREO 
WRITEQUTPUTTAPE NT2, 7, FREQUJFREQ), NSTRIP 


D{ i9 ISTRIP LeNSTRIP 

CGIRTISTRIP, 3) = Of 

CGIRCISTRIP.4) = Of 

ALORR = BLISTRIP) 68 

IF (FREQ( JFREQ)) 17.17.18 

CGIRCISTRIP,1) = ONLOBIOBREDELYCISTRIP)/S 
GIRTLTSTRIP se 2) . O. 

GOTO i’ 


EKJT =1./FREQTJFREQ) eB 1OBR 
CALL RESSEL ( EKJTeLeEJOX-EVOX,0) 


Dt N (EF JOX(Z)¢EVYOR(L) Dee2e¢(EVOR(2)-EJON( 1) Dee? 
f (FE JOKM( 2, © (CE JOK(2) CEVYOR(C LID CEVOXRO 2) eC EVYOR(2)-EJOX01)))/70EN 


G=~ CE VYORC7) ©E YORE LI CE JOK02) eb JOKT1L))/7DEN 
PHIKR = EJOXK(L) oF +E JORI2Z) 06 

PHIKE = EJOK(1L) ©G-E JOX(Z) oF eE SOK?) 

TEMP = ExJlex(ISTRIPIZBIISTRIP) 

EPSR = COSF (TEMP) 

EPSI =-SINF (TEMP) 

CON? = CONL SRL OBReDELYETSTRIPI/S 
CCIRCISTRIP, 1) <CON2 © 1 PHIKROEPSR-PHIKISEPSI) 
CGIRULSTRIP. 2) *CONZ © ( PHIKROEEPS ICP HIKISEPSR) 
CONTINUE 


WRITEOUTPUTTAPE NIT2, 9 


l=LeNSTRIP) 


T/ils 


HMO 200860 
HMOZOU041 


HMO 20046 
HMO200467 


HMO20051 
HMOZO0S2 
Hm™O20053 


mmOz70US4 
4mO20055 
mm020056 
HMOZ005S7 


HMO2Z20058 
HMOZO0O5Y9 


4020061 
m0 20062 
HMOZ0063 
HMO200606 


HMO2Z20065 


4M020069 
HMO2Z0070 
HMO20071 


HBO20075 
HMO20076 
HMOZ00TT 
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STORAGE NOT USED BY PROGRAM 


oec Cl ofc 8=6a FF 
1496 02730 32558 717456 


STORAGE LOCATIONS FOR VARIABLES APPEARING IN COMMON STATEMENTS 


DEC oct DEC ocr | DEC oct orc ocr 
NT1 3756) T?461 NI2 3275@> 77466 =2°3 mT3 32559 77457 P ; 
STORAGE LOCATIONS FOR VARIABLES ABPEARING IN DIMENSION AND EQUIVALENCE STATEMENTS 
DEC Oc! DEC oc} DEC oct DEC OC I 
B 1295 O2417 LGR B95 OlLS7? DELY i595 O25S65 es0n 1095 O7TO? 
Eyox 995 O1743 FREQ 1495 O2727 x 1.195 02253 


STORAGE LOCATIONS FOR VARIABLES NOT APPEARING IN COMMON, DIMENSION, OR EQUIVALENTE STATEMENTS 


De ( OCI DEC Cl DEC OC T Dec OCT 
BIUBK 495 OO7S7 BR 0=s« & 9H -OO TSS CONL #93 OOTSS COnz #92 OO7S4 
COSLAM 491 OO753 DEN 490 00752 eKJT 86694869 GOTSI = = o0750 
ePSe 4A? OOTS? F #86 00746 G #65 OOT45 00744 
TTWO 483 00743 NCARD 482 OCOT4? NFREC 4B1 OOT4I1 NPUNCH 480 OOT40 
NSTRIP 479 OOT3? PHIKI 478 OOT7 Sé PHIKR 477 00735 S 476 O07 34 
TEMP 475 00733 TwOP] 474 O07 32 


SYMBOLS AND LOCATIONS FOR SOURCE PROGRAM FORMAT 


EFN LOC EFN LOC EFN LOC cFN LOC 
8) 1 1 00722 8)2 2 00720 8)4 4 OO716 8)5 5 00657 
8)6 6 00637 8)7 1 00626 8)9 9 00604 B)A 10 00547 


LOCATIONS FOR OTHER SYMBOLS NOT APPEARING IN SOURCE | 


DEC ocr DEC OCT DEC OCT DEC OCT 

1) 467 00723 2) 342 00526 3) 3447 00533 4) 32767 Trill 

6) 350 00536 C)G1 471 OO7T27 C)IG2 472 00730 C)100 473 OO731 
bIz0F 165 oO745 0) 30F 164 00244 0140EF 145 GO221 


c 


LOCATIONS OF NAMES IN TRANSFER VE 
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L1T1TOZ2Z0WH 
STICZOWH 
TTTOZ0WH 
C1LIOZ20WH 


+ 
- 
z 
~ 


SOTOZOWH 
SOTOZ0ne 


*) *x) OATH WIV) 


T's 


{ Ad 


NOT SNIWIO 


‘(lord 


(L)A4 


660020nh 


86002 0WH 


NILNOWANS 


138538 3 


S*f AINO 31N¢@WO) 


» a « a> & my 
zr as i 

o° oo bik 
mm RO La) nN my f 
oO Oo oOo O Bil) 

oo o © 1 ha yh 
oo oo it 
ww ~ ~ oO Oe 


BESSEL 7/09/62 PASE 3 


Y 


1) 


STORAGE NWOT USED BY PROGRAM 


§ -* 


DEC ocr DEC ocr =. 
a25 00505 327561 77461 
STORAGE LOCATIONS FOR VARIABLES NOT APPEARING IN COMMON. DIMENSION, OR EQUIVALENCE STATEMENTS 
DEC ocT DEC oct DEC oct DEC ocT 
324 00504 i 323 00503 
LOCATIUNS FOR OTHER SyMBILS NOT APPEARING IN SOURCE PROGRAS 
DEC oct DEC ocl DEC oct pec acT 
318 o04T6 2) 791 oo84) 3) 795 oOoOeeT 6) 305 99961 
S310 00466 


BJYO 


LOCATIONS OF NAMES IN TRANSFE? VECTOR 


DEC oct DEC oc! DEC oct DEC oct 
9 90000 BJY1 1 90001 


ENTRY POINTS TO SUBROUTINES WOT OUTPUT FROM LIGRARY 


BJY1 
EXTERNAL FORMULA NUMBERS WITH CORRESPONDING INTERNAL FORMELA NUMBERS AND OCTAL LJLATIONS 


I-N OL OC EFN IFN LOC ern 8 6CTPR COLO tru =O 1 FN 
? 00164 12 10 00176 14 11 600262 14 12 55257 
13 00212 8 16 00217 20 15 00226 22 16 90239 
17 00247 26 18 00252 28 19 00271 40 20 00275 
23 00325 4 * to to 273 oOssT 8 oe anaae 
27 00353 “2 30 as 3h 004625 %6 32 00 


00001 
0000 00 QO OQO0000 
“0654 00 & 00000 
0500 00 4 00001 
0671 OO 0 00073 
0500 00 & 00003 
0621 OF 90027 
0500 00 « 00004 
0671 00 O 00030 
0500 00 4 90002 
0621 00 O 00017 
0500 00 0 O 10) 
0601 %0O O C0 
0100 00 O udu 3 
~O0120 Of 93 00071 
0500 00 O 000.6 
0621 00 0 00027 
0020 00 0 00023 
0500 00 O 00046 
0621 00 O 00030 
0500 00 0 00000 
0074 00 4 00053 
6 90000 0 00000 
0020 00 0 00033 
0601 00 O 00000 
-0600 00 O 00000 
~0534 00 & 00000 
00270 00 « GO0OS 
0560 00 0 00037 
0500 00 0 00060 
0020 00 0 OOO2T 
0 00000 O 00041 
ea77TrTiirrririrrT 
#201400000000 


-0634 00 1 OO112 
~0634 00 2 001046 


8JYO 


an | 


B2 


83 
CALL 


B7 
88 
89 
COMMON 


BF 2F 


ENTRY B8JYO 


HTe2 

SxD BIYP-ie* 
CLA 1-4 

Sta 63 

CLA 3,4 

STA BS 

CLA 44 

STA 85 

CLA 264 

STA 61 

CLA 

STO CALLel 
TZt B83 

TMi 82 

CLA B7 

STA BS 

TRA 83 

CLA B?7 

STA B85 

CLA 

TSX BF2ZF.4 
PIE 

TRA 86 

STO 

STO 

LEAD BJYO-1s4 
TRA 5,4 

Loo 68 

CLA B9 

TRA BS 

PZ— VTOMMON 
OCt srvrrrerigtT? 
DEC 1. 

BSS 10 

Sx0 BF2F¢0031,1 
SxD BF2F*¢0025-.2 


STORE ADORESS OF ARGUMENT 
SET J STORAGE 


SET ¥ STORAGE 


SET PARAMETER 
BOTH J AND Y 
OMLY J 
ONLY ¥ 


ERROR RETURN 


BESSEL FUNCTIONS J ZERO AND/OR ¥Y JZERD OF XK 


1/09/62 


Pace i 


4402015) 
44020161 


44920165 
44929165 
44920157 
44020153 
44020167 
ago20179 
4492017) 
4492017? 
44020173 
aeOZ201L7% 
#8020175 
48020175 
44525177 
44920173 
44920173 
4402016) 
4490201451 
ag929182 
44929183 
4402019% 
44920185 
44020155 
#4929187 


44020185 
443279133 


4432917) 
4492911 
@4020192 
#4020193 
aMo720196% 
44929195 


44920195 
449020197 


7/09/62 PASE 2 


SXD BF2F*0050,% 440201935 
i 4492019? 


F2F +00 22 
BF 2F*0009 44920293 
2s% 44992920% 
~OMMON?¢006 440202095 


sans 


LOO BF2F 449202093 
TLQ BF2ZF C0016 44520293 
CLA BF2FeO10?7 449209210 
TRA 434 44929211 
LED SFZF+0050.1 pen a 
BF 2F +0061 1O213 
BF 2F +0042 ei batt 
-OMMON+006 Ft? 44575715 


BF2F e006! 44929215 
-OMMON*000 44920217 
COMMON + 000 5 


0601 


cooe@eeoooreoco-ceeeocooeerrccegecreoeoors+ 


ooTe BF ZF +0088,.4 

6 00000 REDE +0108 4452702721 

0560 OO lel 44929222 

0162 OO BE 2F*eO0029 44929223 
arzFeooT? mmeo2022% 
AFOF e111 4" 227 
COMMON? 009 44927292<3 
~OmMMONF 006 4492022? 

z Z eS ae 

COMMON+ 008 44929236 
COMMON? O009 44920235 


0765 00 O 00043 
0020 00 O OO1T7I 
0500 00 0 00050 
0074 00 4 OOFIL6 
00270 00 O OO176 
0601 00 0 000464 
0500 00 O OO150 
0241 00 O 00050 
~0600 00 0 00041 
0074 00 4 00703 
0 00000 0 00245 
0601 00 O 00046 
0074 OO « 00203 
0 00000 0 00256 
0601 00 0 90052 
0500 00 0 00050 
0302 00 0 00052 
0601 00 0 00045 
0500 00 1 OOOO! 
0340 00 O OO104 
0020 00 0 00164 
#¢202600000000 
0500 00 0 00045 
0074 00 & 00255 
0241 00 0 00066 
0260 00 O 00046 
0560 00 O 00047 
0162 00 0 00160 
0760 00 0 00002 
0601 00 0 00051 
0560 00 1 OOOO! 
0162 00 O OO16% 
0020 00 0 OO172 
0500 00 0 00045 
0074 00 4 00256 
0241 00 0 00046 
0260 00 0 00046 
0765 00 O 00063 


L RS 
TRA 
cua 
Tsk 
TRA 
sToO 
CLA 
FOP 
sr¢ 
Sx 
ele 
$To 
i | 
Pit 
sto 
ae 
Fse 
sTO 
~ i 
AS 
TRa 
pec 
cLa 
TSK 
OP 
+ =p 
Loe 
TOP 
CHS 
STO 
LOY 
TOP 
TRA 
Cia 
TSK 
FOP 
Fup 
LRS 


35 
BF2F¢0078 


COMMON?+OO7 FCs 
BFZFeOL6344 
BF 2Fe0083 
COMMON¢0O03 
BF2F¢0061 
COMMON?O007 
COMMON*+ 000 
BF 2F+0088,4 
BF 2F e022 
-OMMONe 005 
BF 2F¢0088,% 
BF 2FeOL29 
_~URRON+ 009 
COMMON? 007 
COMMON+009 
COMMON? 004 
Ll 

BF 2F¢00?75 
erzFeooTs 
3. 

COMMON +004 FCS 
REIFF *¢O1L 30.4 
~JMMONF0O0 5 
COMMON?¢005 
OMRON? 006 
BF 2F +0069 


COMMON?+ O08 

lel 

BFE 2FeO073 

afzfeoors 

COMMON*006 FCé 
BFZFeOL Sle 4 

“OMMON?C 003 

COMMUN?¢005 

35 


7/09/62 


pace 3 


44920235 
44929237 
449029233 
49020237 
4402024) 
44925261 
44929262 
44020243 
44029245 
ma0202465 
44020265 
449029267 
44929263 
44920247 
m40292509 
49020251 
44920252 
44925253 
4492925% 
44920255 
44070255 
49920257 
44920255 
44952025) 
4492025) 
44929261 
49929252 
44020263 
44020266 
44929265 
44929255 
449020257 
44020253 
qeon02T) 
445275271 
44929272 
44929273 


0500 00 O 00051 
~0534 00 1 OOLIL2 
~0534 00 @ 00155 
~0534 00 2 OO104 

0020 00 4 90003 
-0534 00 1 OOLL2 
-0534 00 2 OO1LOS 
~0534 00 4 OO145 

0500 00 0 00047 

9070 00 4 00007 

0500 00 4 00001 

0621 00 0 002712 
“0754 00 0 OOO0U7 

0534 00 2 00205 

0020 00 0 00212 

0560 00 0 00042 

0260 00 0 00041 

0300 00 2 00000 

0601 00 O 00042 

2 00001 2 002710 

0020 00 *& 00007 
*#105447119564 
+200400000000 


¢164670315620 
+171402400064 
+174554074002 
~ T7303765545 
#201503777346 
202637777765 
+201400000000 
~165404416744 


#171430340621 


200466040575 
*177570222373 
#16445 7452673 


DEC 


DEC 


~OMMON#008 
BF2FeOO31, 1 


3,4 
REZFeOO3lL, | 
BFS ¢0025e2 
Br 2F 4005054 
~OMMON*4 006 
2e% 
i,s¢ 

BFE 2Fe0095 

7 

BF 2*+0090,2 


v2 
COMMON? O01 
BF2Fe¢OO93e2e1 


#.00027100.-.003946464,4.0466679 


~~ 3163866, +1.26596208,~-2.249999T, +1. 


~-00024846,+4¢.00427916,-.0426121% 


+. 25300117. ~~ 74350384, +.60559366,+ 36746651 


#.00014476,-.00072805,¢.00137237 


7/09/62 


PASE 4 


4490 29027% 
44929275 
44920275 
44029277 
44920273 
44929273 
4492023) 
44929231 
44020232 
44929283 
44902028% 
44929235 
449202935 
44929287 
44920283 
44020237 
4992029) 
449259291 
44920292 
44920293 
448920295 
44929295 
44920295 
44925297 


44020293 


44929297 


44920301 


vt 


~166575552336 
#167547604 3046 


~163616754512 
~L?7L552174610 
=1546553077T2 
#2006 30610514 
~1646434251725 
#16546 3450260 
#16645 39617541 
~1L70530121603 
#1625135537003 
174525237266 
#20062 27077324 
#*201622077325 
0300 00 0 OO31% 
0240 00 O OO314 
-~0600 00 0 00061 
0502 00 O OO041 
-0634 O00 & 00063 
“0534 O00 & OO267 
~0300 00 0 OO31L5 
“0773 00 O OOOILD 
0760 00 0 OOO010 
0760 00 O OOOO! 
1 000046 & 00270 
0300 00 O 00315 
0300 00 O 0004) 
Z 000046 & OO274 
0760 00 0 00002 
0601 00 0 000641 
-0300 90 0 00315 
0260 90 0 00061 
0601 00 0 00047 
0560 00 0 00310 
0260 00 0 000462 
0300 00 « OO315 
0765 O00 O 0004635 
2 00001 *& 00301 


Dec 


DEC 
FAD 
FON 
STO 
CLS 
SxO0 
i. XO 
VFA 
ROL 
RND 
LAT 
Tx] 
FAD 
FAD 
Tix 
CHS 
S$To 
Jt a 
PMP 
STO 
Lde 
FmPp 
FAD 
LARS 
Tix 


~.00009512,-.00552760,-.00000077,+.797884656 


~.00013558,¢.00029333,¢.090054125 


~~-00262573,. +. 000039546, ¢.94165597,¢. 78S 39816 


1.57079633 
BF2ZFeO161 

BF 2FeOLSl 
COMMON? 000 
~OMMON#¢000 
~“OMMON*00?,.4 
BFE2FeO140,% 
BFZFe91L62 

La 


BE ZEeCOLSL - Se 
HF ZF eOL62 

L OMRON O00 

BF ZFOO1L 65, 4,4 


~OMMONe O00 
BF2Fe¢O0162 
COMMON?¢000 
COMMON?# 001 
BFZFe¢O15S7 
COMMON? O01 
BF2Fe0O162,4 
35 
BF2Z2F¢0150,4,1 


400 PI/2Z FOR COSIVE 
OBTAIN NO. OF QUADRANTS 


FIX FOR QUADRANT 40 55ST 


OBTAIN WEAREST SEMICIROLE 
TEST QUADRANTS 
FIRST ANI 


PLOAT 


SECOND AND 


JBTAIN SQUARED 


FIJ/aTH+ 


THIaD 
FIRST QUADRANT EQUIVALENT 


ARGUMENT 


7/09/62 


44029392 


449293)3 


44020308 


44920305 
4029395 
44020397 
44929398 
449290359) 
4492931) 
44929311 
*4070312 
44929313 
44029314 
44929315 
449279315 
449290317 
*4029519 
44920319 
44929529 
44929321 
4492932? 
44929323 
aa02032% 
44929325 
44929325 
4492593277 
44929323 
449029327 


00305 02760 00 0 O0O00641 FeP COMMONS O00 
00306 -0534 OO & 00043 LxD COMMON?+002,%4 
00 4 00001 TRA 164 


00312 +¢175506321703 
00313 -200512567414 OCT 600512567414, 201622977324, 234000900000 
00314 *#201622077324 


STO COMMON? O00 . SAVE ARGUMENT SIGN 
BF2FeOL83 TO ExIv IF ZERO 
COMMON O01 Save N 
BF2F¢0186 PUTE TAILAL WALUE, & 

2 x 

A000 COMMON?*O0L x > === 
i x 
REDE eOLAaT 
~OMMON#¢000-% 00 S8ETURN ADDRESS 
BF2ZFe¢O164,% S: INDEX FOR 3 ITERATIONS 


“OMe 
cCOmmONeO0S x N/K 
— pMMONF002 xe 


SKIP ONE 
DO NWOT SKI? ONE 


00343 0162 00 4 00002 
00344 0020 00 « OO00!1 


COMMON+¢000 
Ble cas MASK EXPONENT 


00351 0601 00 0 00061 
00352 -0320 oo 0 00611 


a: 


7/09/62 PASE 6 


4492033) 
449029331 


4492933% 


44 


337 


44929333 
4492033? 


44)23934% 
44920365 


349 
4492935) 
44920351 


* 


44929355 
44920357 


49929351 
44020352 


00415 
00416 
00417 
004270 


0300 00 O 00041 
0765 00 O 00033 
0767 00 0 OO001 
0402 00 0 90412 
-~0501 00 0 00421 
-~0600 00 0 00042 
0300 00 0 00421 
0560 00 0 00042 
0601 00 0 00042 
0500 00 0 00413 
0763 00 0 00033 
0300 00 O 00417 
0601 00 O 00041 
0302 00 0 004620 
0240 00 0 G0041 
“0600 00 0 00041 
0260 00 0 00041 
0601 00 0 00043 
0560 00 O 004156 
0260 00 0 00043 
0300 OO 0 00415 
0765 00 0 000463 
0760 00 0 00043 
0300 00 0 00414 
0765 00 O 00043 
0260 00 UW 00041 
0300 00 O 000462 
0765 00 0 00043 
0260 00 0 00422 
0620 00 4 00002 
¢377000000000 
0000 00 0 004601 
0000 OO O 900700 
#202561251001 
#200754215603 
+200467532521 
+¢20055202 3631 


#201552023631 


~OMMONe 000 
27 

1 . 

BF2F 40223 
BF2F +0230 
COMMON?O001 
BF2F¢0230 
COMMON? 001 
COMMON?001 
BF ZF +0224 

27 

BE2FeO228 
~OMMONe 000 
BEF 2F¢0229 
COMMUN *¢ 000 
COMMON+000 
COMMON +000 
“OMMON* 002 
BF2Fe¢O0227 
~OMMON#002 
af 7% *¢02776 

35 
COMMON+002 
BEIF 40225 

35 
~OMMON#000 
COMMON? O01 
35 

BF 2F +0231 
2.4 
377000000000 
257 

i268 
2.8853912903 
«9614706323 
~59B9TH6496 
- 707106781187 
1.41421356237 


NORMALIZE 


7/09/62 


PASE 


? 


44929353 
4492)935% 
44929365 
44929365 
44920357 
44925363 
44902035? 
4492937) 
44020371 
44020372 
44020373 
4492937% 
44929375 
44929375 
ae07z0377 
44920373 
qagz20377 
44932335) 
449209381 
449203862 
44020383 
4492036% 
44020385 
449293865 
44929387 
44929353 
44920357 
44992939) 
44920391 
44927939? 
44920393 
4492039% 
44020395 
449029395 


44020397 
44929393 
44920397 
449294)) 


00415 
00416 
00417 
004270 


0300 00 O 00041 
0765 00 O 00033 
0767 00 0 OO001 
0402 00 0 90412 
-~0501 00 0 00421 
-~0600 00 0 00042 
0300 00 0 00421 
0560 00 0 00042 
0601 00 0 00042 
0500 00 0 00413 
0763 00 0 00033 
0300 00 O 00417 
0601 00 O 00041 
0302 00 0 004620 
0240 00 0 G0041 
“0600 00 0 00041 
0260 00 0 00041 
0601 00 0 00043 
0560 00 O 004156 
0260 00 0 00043 
0300 OO 0 00415 
0765 00 0 000463 
0760 00 0 00043 
0300 00 0 00414 
0765 00 O 00043 
0260 00 UW 00041 
0300 00 O 000462 
0765 00 0 00043 
0260 00 0 00422 
0620 00 4 00002 
¢377000000000 
0000 00 0 004601 
0000 OO O 900700 
#202561251001 
#200754215603 
+200467532521 
+¢20055202 3631 


#201552023631 


~OMMONe 000 
27 

1 . 

BF2F 40223 
BF2F +0230 
COMMON?O001 
BF2F¢0230 
COMMON? 001 
COMMON?001 
BF ZF +0224 

27 

BE2FeO228 
~OMMONe 000 
BEF 2F¢0229 
COMMUN *¢ 000 
COMMON+000 
COMMON +000 
“OMMON* 002 
BF2Fe¢O0227 
~OMMON#002 
af 7% *¢02776 

35 
COMMON+002 
BEIF 40225 

35 
~OMMON#000 
COMMON? O01 
35 

BF 2F +0231 
2.4 
377000000000 
257 

i268 
2.8853912903 
«9614706323 
~59B9TH6496 
- 707106781187 
1.41421356237 


NORMALIZE 


7/09/62 


PASE 


? 


44929353 
4492)935% 
44929365 
44929365 
44920357 
44925363 
44902035? 
4492937) 
44020371 
44020372 
44020373 
4492937% 
44929375 
44929375 
ae07z0377 
44920373 
qagz20377 
44932335) 
449209381 
449203862 
44020383 
4492036% 
44020385 
449293865 
44929387 
44929353 
44920357 
44992939) 
44920391 
44927939? 
44920393 
4492039% 
44020395 
449029395 


44020397 
44929393 
44920397 
449294)) 


POST PROCESSOR ASSEMBLY DATA 


#23 1S THE FIRST LOCATION NOT USED BY THIS PROGRAM 


REPERENCES TO DEFINED SYMBOLS 


442 R 
l2 Bl 
21 B2 
23 B3 
a7 84 
30 B5 
33 Bé 
36 BT 
37 88 
40 Bg 
53 BF 276 

1 BJYO 
2 8 86—C ALL 
4! COMMON 


11 
15 


14, 
i?, 
22 


21 


53, 
104, 
137, 
204, 
320, 
£02, 

31 


64, 
136, 
210, 
324, 
367, 


NO ERROR IN ABOVE ASSEMBLY. 


20 
35 


54, 
106, 
140, 
206, 
3226 
407 


66, 
L4l, 
2ll, 
327, 
371, 


$5, 
LOT, 
146, 
207, 
326, 


76. 
213 
S31, 
372, 


ST, 
Lil. 
147, 
214, 
$330, 


100, 
143, 
257, 
332, 
373, 


él, 
Li5. 
156, 
256, 
340, 


102, 
151, 
261, 
334, 
376, 


62, 
116, 
162, 
262, 
352, 


110, 
153, 
271, 
335, 
401, 


67. 
Liles 
163, 
263, 
456, 


113, 
154, 
274, 
336, 
404, 


Tle 
124, 
172, 
270- 
361, 


121, 
160, 
27T, 
342, 


7/09/62 


73, 
126, 
173, 
272, 
36%, 


122, 
L64, 
s01l, 
351, 


Pa 


SE 1 


7S. 
ti, 
175, 
eer 
37D, 


135, 
LoT, 
$26, 
360, 


7? 


134 
177 


302 
375 


132 
171 
316 
362 


> > 


00041 


00053 


0000 00 

-0634 00 
6500 00 

0621 00 

0500 00 

0621 00 

0500 00 

0621 00 

0500 00 

0621 00 

0500 00 

0601 00 

9100 00 

“0120 00 
0500 00 

00 

00 


@wmecoccoeFoococo#ogocco9eoeqoccoeo# or or 9 


-~0634 00 1 


00001 


00027 
00004 
00030 
00002 
00012 
00000 
00075 


BJY1 


61 


B2 


CALL 


COMMON 


BF SF 


ENTRY BJYI1 


HTR 
Sac 
CLA 
STA 
CLA 
STA 
CLA 
STA 
CLA 
STA 
LA 
STO 
Tze 
Tel 
CLAa 
STA 
TRA 
CLA 
Sta 
Cia 
TSX 
PIE 
TRA 
sTo 
ste 
Lxo 
tra 
LB 
CLA 
TRA 
ere 
oct 
Pre 
BSS 


$xD 


BESSEL FUNCTIONS Js 


BJYin-ie® 
1,4 

83 

44% 

B4 

&,* 

65 

2,54 

Al 


CALL el 
6% 
B2 
87 
AS 
83 
a7 
6» 


oF aF.* 


86 


BJYi-1+4 
5.4 

88 

a9 

ws 
-ORMUN 


STTTTTIITIATT 


° 
10 


7/09/62 


(K) AND Yo (8) 


STORE ADDRESS 
SET J STORAGE 


SET ¥ STORAGE 


SET PARAMETER 
BOTH J AND Y 
ONLY J 
ONLY ¥ 


FRE OR RETURN 


JF ARGUMENT 


BESSEL FUNCTIONS J ONE AND/OR ¥Y INE OF XK 


BF 3F*¢0033,1 


PASE 1 


449206095 
449206)7 
44990206098 
449020460? 
4402061) 
44929611 
44929612 
44)2)413 
440Z0615% 
440279815 
49929615 
44929617 
44020613 
44920617 


449529423 
4492062% 
44929625 


440209425 
449279462) 
440270643) 
44920631 
aeoz0e3? 
229433 
oe3s6 


4452734635 
44029635 
44920637 


44975661 
440294642 
449294643 


0260 00 
0601 00 
0974 OO 
0 900000 
9765 OO 
0760 OO 
0560 00 
O162 00 
60620 00 
0601 00 
0074 OO 
9 9oo00 
O241 00 
~0600 00 
5560 O00 
0260 00 
00746 OO 


rcoceeercoewaocoroeoooooreccoe@e-cceeer~ 


00047 
00153 


TRA 
$10 
TSK 
Pze 


stg 
LDQ 
FMP 


TSK 


BF3F*¢0025,2 
BF3Fe0053,4 
Lee 


BF 3F*eO103 
BF 3F +0009 
2.4 

COMRON+ 006 


COMMON+0O07 
BF AF eO1 03 
BF SF eOO16 


34% 
6F 3 +0053,1 
BF 3F +0064 
BF $F *00645 
~ OMMRON* 006 FL2 
BF 3F e006% 
ne 
OmrOnN*+000 
COMPON +000 
BFF +e0097,4 
BESFeO1 11 
45 
COmRRON+ O06 
1-1 
BF $F e003! 
BFF e008) 
Cc OMRON O08 
BF 3FeO0092,% rea 
BF sFeOL 16 
COMMON? 006 


= 


COMMON +009 
~OMMON+006 
BF3Fe¢O110 

BFSFeOL9L.S 


44) 72 )04%% 
440 204%%) 
aa92046%5 
44020667 
449294%3 
4492944? 
4492045) 
49920451 
4490729457 
449020653 
49929454 
44529655 
44929655 
449279657 
74979855 
44992065? 
449027945) 
44520651 
44929452 
44929653 
q4092085% 
49920465 
489209465 
44970657 
44929653 
4492066? 
aso708e79 
44997294671 
44929672 
445704673 
4492967% 
449929675 
4a4a9298TS 
a4020677 
44029673 
4492729677 
4492943) 
449296351 


001242 
00123 
ooLze 
00125 
00126 
00127 
00130 
00131 
00132 
00135 
00134 
00135 
00136 
Oo137 
00140 
oo14l 
oo14se 
00143 
001464 
001465 
oo14s6 
00147 
ooLso 
00151 
00152 
00153 
00154 
00155 
OO15S6 
00157 
00160 
OOlel 
oole? 
co0163 
OOLes 
00165 
00166 
ooie?7 


0601 90 
0500 00 
0241 00 
-0600 00 
0074 00 
0 90000 
0601 00 
0074 00 
0 00000 
0601 OO 
0500 00 
0302 00 
0601 00 
0500 00 
0340 00 
0020 OO 


00001 
00104 


00167 
+202600000000 


cow oaecoocoroogrcoe@gerocedagco 


00045 
00261 
000464 
00046 
00047 
00163 
oo002 
00051 
00001 
00167 
O01 76 
00045 


oooreocgcoccoCgce# eS 


BF 3FeO00R7 
BF 3FeOL 32 
COMMON+008 
COMMON? 009 
35 

BE SF eO082 
COMMON?#007 
BF 3Fe¢O1L66-% 
BF 3FeO087 
COMMON+00 4 
BF 3F 00646 
~OMMON*007 
~ OMMON#000 
BF 3Fe¢009244 
OF SF 90129 

- 

OMMON+005 
BF SF *0092,4 
BF 3FeO1L 32 
COMMON#¢009 
~OMMONC O07 
COMMON? 00% 
COMMON+004 
el 
OF3F+¢O0025 
RE AF eOOTS 
}. 

COMMON+ O04 
BF3FeOL34.% 
COMMON+003 
COMMON? O05 
COMMON?¢006 
AF 3F e007? 


COMMON+O008 
isl 
REE *eOOTS 


AF 3F e008 
COMMON? 006 


FCS 


FLO 
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44929682 
44929683 
490729088 


445259647 
449296433 
44920657 


3 


72 
44020673 
44020675 
44929695 


44979699 
449209509) 
44020501 


44929595 
499295097 


Pil 
qt waitin 


rt 


' 


os = -_- = : = 


; I 
44929512 
44920513 


44929513 
44929517 


2S = = — 


[> 


00170 0074 OO 4 00260 
ool/7t 0241 00 0 00044 
00172 0260 00 0 90046 
OO1L?3 OF760 00 0 00002 
OO174 0765 00 0 00043 
00175 0500 00 0 00051 
O0O1L?é -0534 OO 1 OOLILS 
Ooo1L7?7 -0534 OO & 00140 
00700 ~0534 00 2 00104 
00201 0020 00 4 00003 
00202 ~0534 00 1 OO114 
Ovi.vus ~0534 O00 2 OO1L04S 
00204 -0534 O00 4 OO140 
00205 0500 00 0 OOO047 
00206 0020 00 4 00002 
O0O207 0500 00 4 O0001 
00210 0621 00 0 00216 
00211 -0754 00 O 00007 
00712 0534 00 2 OO11 
00213 0020 00 O O0OO716 
00214 0560 00 0 00042 
OO715 06760 00 0 OOu41 
00216 0300 00 2 00000 
00217 06601 00 0 0064? 
00240 2 ooool 2 002714 
00221 0020 00 46 30002 
00222 #105447113564 
00223 *16056407457T7 
00274 ~16551502 3444 
00235 +171442421130 
00246 ~174503736715 
OO227 #1766577 761L11 
00230 ~2008377T7TTTS3 
00231 +200400000000 
oozs2 #170555254422 
~174510365414 
*LTT477T7TLL001 


“201521012012 


BF 3Fe¢O1L33,% 
~OMMON#003 
~OMMON+005 


35 
COMMON+008 
AF 3Fe0033,1 
BF 3F*+0053,.% 
BFS +0075 62 
344 
BF3F+OO3351 
RF 3F+0025,2 
BF 3F*¢0053,4 
COMMON+006 
Z+* 

i,* 

BF 3F +0099 

7 

BF 3F e0096,2 
AF 3F 40099 


COMMON+ 001 

BF 3FeOO9T. 251 

2e% 

it-is 

+. 00001 109,~—. 00031761 «*~004643319,-.03954289 


¢.2109 459 73--.-56249985 


> 
¢.0027873,-.0400976, +. 3123951,-1. 3166827 
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4492952) 
449290521 
440720522 
44929523 
449020525 
44325525 
44920525 
44929527 
44970574 
4992052? 
449020535) 
449529531 
4492953? 
44929533 
44070555 
44929535 
44979535 
44929537 
44920534 
#4920539 
#907954) 
449929561 
449029056? 
44929543 
14929548 
44929545 
44020565 
44979567 


44929543 


4992956) 
4402055) 


00236 *#202425423632 
00237 ¢176705022433 
00240 -200505746037 
00241 


-L70507023776 
164546557423 
*LT3I41LT726752 


~ 166642403726 
- 166605630311 
171642024531 
~1L62731 5 
“195 767 
+200622077324 
201622077325 
0300 00 0 00317 
0240 00 
~-0600 
0502 
-0634 
-0534 
0300 
0773 


07T60 


t 
6360 606 


0300 OO 
2 00006 


ccooe~eercoe~eecocor#ogoc°o 


= — [SS ——_ 


= = a yroerve2 pase 5 


DEC #2.1682709,¢.2212091.-.636%198 49929551 


49029552 

DEC #.01659667,%. 00000156, ¢. 79788456 44029553 

DEC = 9070554 

DEC --00005650,.~.127499612+*. TES59816 = ANOTOD>> 
—— = ee 

o€C 1.57579633 | 4E575S55 

FAD BF3FeO164 ADD PI1/2 FOR COSINE 44020557 


BF SFeO1L6S OBTAIN NO. OF QUADRANTS 44029555 


— —— = = =e at — 


=* - = = —=_— : == 


) * 
BFUFe¢O143,4 44373552 
BFE 3F*eO0165 FIX FOR QUADRANT AD UUST 44)29563 
4 4492956% 


4402056? 
4402957) 


Common? 000 
AF SF eCO1L48. oe 


— a == = 


- 


a 


STO COmMMON?O0!1 44920575 


"a 
_* 


0260 00 0 00042 
0300 00 @& 00320 
0765 00 0 000463 
2 00001 4 00304 
0260 00 0 00061 
“0534 00 4 00043 
0020 00 & 00001 
164476053726 
~1714622466361 
*#L75506321703 
~20051 2567414 
201622077324 
#2 34000000000 
0601 00 O OODS6!1 
0760 00 0 00003 
0100 00 0 00445 
0601 00 0 000462 
“0320 00 O OO ¥50 
0765 OO O 00001 
0400 00 O 00042 
0765 00 0 00001 
0400 00 O OO*S1 
—~0634 00 & 00041 
0534 00 4 00322 
060! OO O 00063 
0500 00 O 00042 
0240 00 0 90043 
-0600 00 0 000646 
0500 00 0 00044 
0300 00 0 00043 
0402 00 O 00350 
2 OOOO] & OO334 
“0534 00 46 GOUOSI 
0560 00 0 00041 
0162 00 4 00007 
006020 50 4 OOOO! 
001000000000 


100400000000 


FMP (_OMMON?O0]1 
FAD BF3F *e0165,4 
LaS 35 : 
rae fomonreee 
LxO 2OmMON?+002,¢4 
TRA 1,4 Exit 
OCT 164476053726, 571462246 361,175506321703 
, 
OCT 600512567414, 201622077324. 234000000000 : 


La0 


Log 
tor 
TRA 


DEC 


~OmMRONe 000 SAVE A&GUMEMT SIGS 


" 
BF SF rOLeS TO Exit tf Zeno 
COMMON+001 SAVE © 


BF SF eOl ad PUTE FRIAL VALUE, X& 

l x 

COMMON?+O01 x 

I x 

ef 3F*+0190 

COMMON+ 000.4 OO RETURN ADDRESS 
BFUFCOLAT.S SEY INDEX FOR 3 YIVTERATIOINS 
COMMON? 002 Save x 

~OmmONe OO! COMPUTE SQUARE ROOT 
of x Ms 

-OmMMON*003 x 

COMMON +003 x N/K 

“OMMON+ O02 ek 

AF 3FeOL89 IVIDE BY 2 
BF3FeOL77.4,1 REPEAT OOP 


OOORE EXIT ADDRKE.S 


134217728,8657043456 Oe27 EXP -127, 1/72 © ERP -64 


49)20577 
4499020573 
49920577 
4402058) 
499020581 
49920582 
449929583 
44920588 


49929585 


49020585 


449720582 
4452059) 
44920591 
49929592 

“s0720593 

490720594 

44920595 
44979595 
44)20597 
‘99290593 


029597 
4992069) 


++ 
499205597 
499296)3 


449020696 
442920605 
44929697 
49975553 
4492950? 


00415 
00416 0000 00 
0041? #202561251001 


ecoB FP ecoccoeeeoccoeeeoco@eegocoeeeco°eeoors 


TZt 1s ERROR RETURN 
Tm] Lle* > 
$To Cc 


Fad © 
LRS 27 

ALS 1 

SUB BF 3F 40226 


Q COmMmMON?eOO1 
») -omMONe OO] 
A BE BFeO2 27 


SR BFE 3F e034? 
Ow °° OMMONFOD) 
TO COmMmRONe O00 


FMP "“OMMON+ OOD? 
FAD BF 3Fe0229 
Las 35 


FAD 8F3F*0228 
35 
COMMON? O00 
COMMON? O01 
| 


4HTR 128 


DEC 2.88653912903 


———— — 
—- _ 


= - — 
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4992)951 
4020511 
44920612 
#4929513 
449929515 
449290515 
44D29615 
449929517 


44929621 
4492962? 
44920523 
4070675 


14575577 
44929623 
44)790677 


errr 
49929533 
44920538 
449020535 


49 
4402964) 
44929661 


449029665 
440206% 


= -~ 


00420 #200754213603 
00421 #2004672532521 


0042726 #20054271027T7 


DEL 
DEC 


-9614710632 3 
5989786496 
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449020663 
4492064) 
4492065) 


44525653 
44902065% 


9F 


wou ON 


“A VGW3ISSY JA08V NI 


"oa 
*stl 


‘ce *20€ “TOE “412 8422 «8492 «E92 «8292 «8412 «8612 «HITZ «#502 
‘E91 “COT “261 “961 88ST) fLOl) 89M) 868 foet fTeT 8 9E 
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Cle 
tit 


“CO “ELE “TLE %49€ $96 *Z79€ ‘*19€ “S56 £9 +298 "tte “1k “Sze “Ect 


2c$ 


4 


<— 


; 


ha 
i} 


came wa 


yl vt 
1) EE 
i Ha 


vit 
bie tit 
Hop Pai oat 


Hit 


w¥¥9ONd SIMA A@® O3SN ION NOTLVIO) ISWI14 3HI ST 92% 


Mh 
We 


£y 


5 


6 


1b 
‘ 


SUBROUTINE AEROPS (VBRW,ISTRIP, CH, NTAPES,NITAPET,NSTART 


DI MENS LON CH(100,4), Al22) 


FORMAT {¢ LHO 40K, 24H it«, 
10H THRU HMO? 114 ) 

FORMAT ( 214, 64K, S4HHMO?2 L146 ) 

BE DZ *60 3044000260 : 

1ORG-1 

TS=NSTART 

WRITE OUTPUT TAPE NIAPE?7, 1, VBRw, IS 

[S=1iSel 


K=Ze 1 SIRIP 

WRLTE OUTPUT TAPE HTAPET, 3 Ke ISTREPs IS 

ISs1Sel 

D6 & t#1,22 

A(1)=0. 

Le=C 

M=Q( 

DO 12 Iie LSTRIP 

BO il seis? 
Meme, 
IF (M-23) 6.5.5 
MeM-22 
CALL BINPU (Ae22-10RG, BCOZ-ISeNTAPET) 
LORG= LORG+ 272 
LS=1iSel 
po 6 WNe1,.22 
AIN)=0. 

LF (m-23) 15-55 

if {t) 8,8.,7 
A(M)=CNtI,J) 
MeMel 
If (™-23) 10,99 
L=1 
GOTO 5 

i=0 


) 


1/09/ 


44020655 
4HMO20656 


4M029661 
HMO 20662 


48029666 
44020667 
HMO 20668 
4MO020669 
‘HMO20670 


4mO70671 
HMO25672 


HmMO02067 35 
HMO2067% 
"mB 700 /* 
HMO20676 


wm020677 
HMO2D4678 


48020679 
44020680 


440290685 
HM9 20686 


HHO020691 
HMO2Z20692 


= => od .- 


er - 


NdNI@ WIV) 41 
1¢St=S1 
2-New ee 


(Zor *TIHI=(wdY 


(L3dWIN°S1*7098'0*0'W! 


i iB 
eel Ah 
} i ‘hi iH 
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an 
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| 
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i 
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COLO20WH 
669020WH 
66 9CZ20WH 
€69020WH 
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STORAGE NOT USED BY PROGZAM 


oec 8 6c ogc 86 «Oc 
269 004615 32561 T7461 


STORAGE LOCATIONS FOR VARIABLES APPEARING IN DIMENSION AND EQUIVALENCE STATEMENTS 


DEC oct DEC ocr bec Oct oeCt oct 
a 268 00414 


STORAGE LOCATIONS FOR VARIABLES NOT APPEARING IN COMMON, DIMENSION, OR FOUTVALENTE STATEMENTS 


DEC OCT DEC OCT DEC Oct DEC Jct 
BCDZ 246 00366 LORG 245 00365 is 244 00364 ® 243 90363 
t 242 00362 S 261 00361 


SYMBOLS AND LOCATIONS FOR SOURTE PROGRAM FORMAT STATEMENTS 


EN LOC EFN LOC EFN LOC EFN LO. 
$21 1 90354 8)2 2 00347 6)3 3 OO327 


LOCATIONS FOR OTHER SYMBOLS NOT APPEARING IN SOURCE PROGRAM 


DEC Oc T DEC OCT DEC OCT OEC Des F 
2) 199 00307 3) 204 0031% 6) 206 00316 9) 237 09355 
CIG2 239 00357 C3100 2740 00360 


LOCATIONS JF NAMES IN TRANSFE2 VECTOR 


DEC OoctT DEC OCT DEC oct DEC ICT 
BINPU 2 00002 (FIL) 1 00001 (STH) 0 00000 


ENTRY POINTS TO SUBROUTINES NOT OUTPUT FIDM LIBRARY 
BINPU (Fis) (STH) 


~~ 
— 
‘?. 


4 
H Mi ii i 
a 


ul i 
i Ai i Hit ht 


i Hi si} 
oe 


ad 14: ud 


z1 


$2200 24% 


tt 


S*Z00 S* 


+! 


19200 8% 


SKOTIVZC) WA20 ORV S¥3GWAN VINWYC4 TWNYSINI ONTONOdS3YYO) HIIM SY3QWAN VIOWHO4 TWNW31X3 


1S 


att 

Hi | 
mits it } 
ee te ; 
‘ Mateal | ie 
Haat hh 

' 
ian 
Pil 
ut 
HN 
; WHT i 

in} 

1h 


Pe ae 
‘ TOL 
yar ih 
HMMS 


! 
i 


ty 
tnt 
: 


HN*n*o*o*o* 


—@1L0Z0WH 


*E€IGVIN JdVA ANPING JATYR S 


t 


L1LO020WH 


‘ 7 
hh 
’ 7 


1 8(gte) 


> 


‘Z21LC20WH 
LILOZOWH 


0109 — 


*Z3dV1IN 3dV1i ANGNT OVIU 


1 


yRe 


SOLOZ0WH 
601020WH 


N1Gw INTINOWENS 


*Z3dVIN) 


*EIAVIN 


ae | 


ines) 


330 


130 


(MAY) 
330 
330 


130 


1oooo0 | 
130 


“ 

_ 

5» co 

x 

a 

> 

- rm 

vw 

= z 

= | oO 

~ _- 

as 

] o ia 

m Leal > Baal 

a) oO mn © 
~ H 

ee 

°o < 
oo Oo . 

on a. v 

o- - 2 

~N 7 © 

> 

» 

i” 


(hSi) 


330 


330 
SINIW3LVIS LYWHOS WYYOONd JDBUNOS BWO4 SNOTIVIOT CNV SWORWAS 


130 


C1 T1200 CAV S¥3GKON VIAWECS TWNYSINI SNIONOdS3¥BOD HAIM S¥IGWAN VINWYOS TYNYI1X3 
o00ce © 
120 


- 
- 


tL REE 
AES Eh 
at 


: 
‘ ¢ 
fi 


SACTIVYV 


ts 


00006 


FRANGFER VECTOR 
60000 743146623460 
00001 746651623460 
00002 745123303460 


323 


00006 0634 00 1 001462 
00007 0634 00 2 00143 


00013 0621 00 006 
00014 -0500 60 4 00002 
00016 0602 00 0 77776 


— == 


BimPU ROUTINE TO WRITE COL BIN CARDS ON 


FIetl 


Pace | 
CO CSS SEE SE SES EEEEEESEEEE SESE SEES EH EEE SESESEHESESEESSESESESESEEESeSeHeEEeHEDZDI CC! 
. 44020723 
. ew wmoz0727 
. Se —— 3? 
* ¥ TO BE PUNCHED) >. Sa 1 
. (NO. wORDS TO PUNCH) 5732 
° (CARD ORIGIN FOR LST CARD) 44020733 
. (SEQ NO. OF LST CARD) 4402073% 
. DECK, LS! AND 2ND GHARACT %> 
7 SSS = =-=-= ss or 
eoemeCONTRARY BE | 
. 4 5733 
© LTEMS MARKED (e) MAY BE DELETED. B8CO ID wILt BE 4402073? 
° UNCHANGED AND SEQ. NOS. WILL BE CONTINUOUS STARTING 4402076) 
. FROM O00. ALSO ORDER WAY SE SwIICHED. | 
. = = = ——__ = 
* THIS VERSION PUNCHES 01 ONLY. = 
. 
PT PTTTTTTTTITITTTITTTTITITTT TT ee ee, 
ENTRY BINPU 440209745 
. == = wmozor.7 
(15S) 
(wRS) 
= = — = 
BINPU SKA 449297465 
SKA 44020746? 
cae 
== sTD === 
= cia LOC OF ARRAY === = 
STA 
CAL® #ORD COUNT 4402075% 
Si* 


END=0 IF TRANSFER CARD 44020755 


ii 


00016 0402 00 0 00325 
00017 0622 00 o C0066 


449209755 


44020757 
00023 0120 00 0 00025 
00024 -0501 00 0 00266 ADD RELATIVE BIT 44029762 
00025 -0501 00 0 00334 7-9, WORD COUNT Ze 44929763 
— a 2enee eee . ‘ee seeeecene ide 

o 
. AND SET CELLS FROM CALLING SEQUENCE. 49929753 
pousessoqucneonsssnnnessaqqneesosoossoscscouecoocoosocosesesesSsSeeee eee. 
00033 0322 00 0 00307 ERA MSKTSK 4402077% 
00034 -0100 00 0 00054 TNZ G2 nO MORE vSk cES 44029775 
. - oe : ‘=? 
00041 -0100 00 0 00043 TnNZ e+2 1S SEQ ND NON-ZERO. 44920786) 
T 0008? -0754 00 0 00000 PD “0 44020781 


a / ~ BO = 4 slele‘ale ' Pace 
00047 0602 00 0 00267 Su SEQNO Save 44029785 
00050 1 777?TTl 4 00053 TxI G5e4-1 440207867 


00055 -~0520 00 O 7T7TTT6 nat ENO 1S wORD COUNT ZERO H4920792 
00056 osee oo 0 00152 TRA Taco MUST BE A&A TRANSFER CARD 44020793 


SS 


BEMPU ROUTLNE TO WRITE COL BIN CARDS ON TAPE. FIBII 


0774 00 
“0754 OO 
Oo7T74 00 
0560 00 
~0600 00 
0361 00 
1 ooO001 
3 00000 
2 00001 
0634 00 
@361 00 
0602 00 


~0500 00 
0560 00 
~0765 OO 
-0500 00 
“0763 OO 
-~0600 00 
~0130 00 
0634 00 
0774 00 
0774 OO 
OTTs OO 
-O7TS4 00 
~O114 O1 
-2 00001 
0767 OO 
0020 00 


OOFNrtFNN FF OK 


CoOorocOorrner Cococoe 


fi 


a 


7/09/62 


PASE 3 


44929794 


ROR E REE ESEOOEOEEEHESE EE EESE EERE EE EEE E SESE SEE ES EEE EEE Seeeeeeeeeeeeee44)2D7T75 


BUILD THE CARD IMAGE. 


q807207T95 


PPPUTTTITITITIT TTT TITEL I TTT T TTT 


NEXT 


COUNT 
ARRAY 


LOCN 


EDLT 


44020793 
AXT 2292 4402579? 
PXD CLEAR AC FOR CHECKSUM. 4492089) 
AXT **,% 449290801 
LBe *e,4 MOVE ARRAY INTO CORE. 44929392 
s1¥Q CEMAGE+24,? 44029893 
ACL CIMAGE*+24,2 ACCUMULATE CHECKSUM FOR BODY. 4902089% 
Tx! ®+1,.4,1 44920895 
TXH OUT sSs** FIMISH WHEN SPECIFIED 44920805 
TIX ARRAY s,2,1 BY NO. WIRDS DESIRED. ( 2,4) 449020807 
SKA LOUNT,«& SEY COUNT FOR NEKT LOOP. 49020803 
ACL CIMAGE ADD IN CONTROL WORD. 49920897 
Sis CIMAGE+! PUT CHECKSUM IM IMAGE. 44020819 

449298511 

OOPS COSHH OOEESEEEE ESSE SESE ESE SEE ESEHESEOS SEES SEH ELESesessescccooeseed4)2 0512 
EDIT THE ILOENTIFICATION FIELD. 44920813 

ROR REESE EEE EES EERE EERE EE EESREEEE SSE RESER COS OR ERE R SE HOH OS eR eee eR Eee TIOZOELS 

44020815 
CAL SEQND 44020815 
LDQ Lt1) 449270817 
LGR Lt 449020813 
CAL BCDIO 4492951? 
LGL 6 #9072082) 
STO IOLCOD SAVE FOR FINISHING. 44929621 
KCL 44020822 
SKA Svi,l 44929823 
AX] Sec 4402082% 
AXxT 20% 44920825 
Axi 3,1 44020825 
PxXD 44929627 
CAQ TABs st 440209823 
TNX e+3,1,1 44920827 
ALS 12 4402083) 
TRA e-3 44929831 


9S 


0602 00 
lL 777TT 
zZ 90001 
0560 300 
3 90000 
0774 00 
0761 OO 
-0500 00 
0074 OO 
0522 60 
~O774 00 
0522 60 
0754 OO 
0621 60 
0074 00 
=~0500 00 
0400 00 
OL1L4 06 
0602 00 
0520 00 
9020 00 
-0500 00 
0602 60 
0774 00 
0774 OO 
0774 OO 
9020 00 


ernerOocccoccoocc#oro# oOFeoor 


BIN AT LAST TO LAST#3 44020832 


490206 33 


Svi 
* 44920833 
SPODereosoecoaseoseoesosessocsoosocosoossossooossoooossossousoooeelons seh?) 
. Fem Is ‘BUTLT: WITH THe BODY a = 
. “AND ID AT LAST THRU LAST#3000 2 2 2 
= = SS 
a PENDS CECT CCECCUCEECCCOORSTOCCSSOSCCCC SC COC CCC OOOSDCOCSEEEOS 
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